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INTRODUCTION 

This article will survey some of the research in the area of 
computational unsteady transonic flows about airfoils and wings, 
including aeroelastic effects. in the last decade, there have 
been extensive developments in computational methods in response to 
the need for computer codes with which to study fundamental aerody- 
namic and aeroelastic problems in the critical transonic regime. For 
example, large conmercial aircraft cruise most effectively in the 
transonic flight regime and computational fluid dynamics (CFD) pro- 
vides a new tool, which can be used in combination with test facili- 
ties to reduce the costs, time, and risks of aircraft development. 

One of the major uses of unsteady transonic aerodynamics is in 
the flutter analysis of supercritical wings. Experiments have shown 
that dips occur at transonic Mach numbers in the flutter boundaries 
for wings and such dips are especially severe for supercritical 
wings. This phenomenon is attributable to the motion of shock waves 
on the wings. The proper modeling of the physics of such moving shock 
waves requires that the CFD methods solve nonlinear partial differen- 
tial equations for regions of mixed subsonic and supersonic flow. 
Currently, the most advanced codes use potential equations for model- 
ing the flow; such codes are being used for generic research in aero- 
elasticity. More advanced codes are under development that use the 
Euler and Navier-Stokes equations; such codes can model vortices in 
addition to the shock waves. 

In this article, first a result ^ 2 ^ will be shown of transonic 
flow with moving shock waves over a pitching airfoil. Then results 
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will be presented for flow over (1) a transport wing,^ including 
aeroelastic calculations and (2) a low-aspect-ratio wing.^1 j n a n 
cases, comparisons with experimental data will be made. 

TRANSONIC FLOW OVER AIRFOILS 

The development of algorithms and codes for simulating unsteady 
transonic flows by solving the unsteady transonic small disturbance 
potential equation was started at NASA Ames Research Center by 
Ballhaus and Lomax. Initially a code, called LTRAN2, was developed ^ 
for a low frequency approximation to the governing equation. Later 
the code was generalized to account for all frequencies and was called 
ATRAN2 . The algorithm used in the code is an alternating-direction, 
implicit (ADI), finite-difference scheme. It is conservative and 
time-accurate so that shock wave motions are modeled correctly. There 
have been many improvements to the code, including improvements in: 

(1) accuracy, such as the use of second-order differencing; (2) sta- 
bility, such as the use of monotone schemes; and (3) capability, such 
as the inclusion of viscous modeling and the capability to model wind 
tunnel walls and supersonic free streams. See refs. 1-2 for details. 

In Fig. 1, a calculation^ is shown for the ATRAN2 code. The 
results are for transonic flow over a pitching NACA-0012 airfoil. The 
calculations show the effects of including viscous modeling and are 
compared with experimental data. Figure 2 shows results over a cycle 
of motion, during which the shock wave appears, grows quite strong, 
and then disappears. This case demonstrates the necessity of solving 
a nonlinear governing equation in a time accurate manner in order to 
correctly account for the shock wave behavior. 

TRANSONIC FLOW OVER WINGS 

The two dimensional algorithm has been extended to three dimen- 
sions and there have been extensive code developments and computations 
of transonic flow over wings. 

The first case shown here is for flow over a transport wing.^ 
Figure 3 shows the model tested in Japan. Figure 4 shows results for 
steady flow (using the code XTRAN3S), which show the presence of a 
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strong shock wave. Figure 5 shows results for unsteady aerodynamics 
with comparisons that show the inadequacy of using linear equation 
methods (doublet lattice) at transonic speeds. Finally, Fig. 6 shows 
results for the flutter boundary in the transonic-dip, Mach-number 
range. 

The second case shows results for flow over the F-5 wing. For 
this wing, algorithm developments were needed for an adequate computa- 
tional grid.^l Fig. 7 shows the model and Fig. 8 shows the pitching 
motion. Finally, Fig. 9 shows the comparison of the computed results 
(using the code ATRAN3S), with experimental data of the unsteady pres- 
sure coefficients on the upper surface. The agreement with experiment 
in both the real and the imaginary components near the midchord region 
of the wing shows that the motion of the shock wave, which occurs 
there, is correctly simulated. 
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Fig. 1 Unsteady pressure coefficients for a NACA 0012 airfoil, 

M = 0.755, k = 0.1628, mean angle is 0.016 deg, pitching angle is 
2.51 deg. 
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Fig. 2 Time history of upper surface pressure coefficients, 
NACA 0012, M = 0.755, k = 0.1628. 
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FLUTTER ANALYSIS OF JAPANESE TRANSPORT WING 



Fig. 3 Code correlation model of 
Japanese transport wing. 
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Fig. 4 Steady state initial 
conditions for the transport 
wing, M = 0.85. 
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Fig. 5 Modal aerodynamics for the transport wing, M = 0.85, 
k = 0.121. 



Fig. 6 Flutter boundaries, computations and experiment. 
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MOOAL MOTION IN EXPERIMENT 
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F5 WING PLANFORM 
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Fig. 8 Unsteady modal motion 


Fig. 7 Dimensions of the F-5 
wing. 


of the F-5 Wing. 
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Fig. 9 Comparisons of unsteady pressure coefficients, computations 
and experiment. 
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